Skyrmions in Spinor Bose-Einstein Condensates 
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We show that spinor Bose-Einstein condensates not only have line-like vortex excitations, but in 
general also allow for point-like topological excitations, i.e., skyrmions. We discuss the static and 
dynamic properties of these skyrmions for spin-1/2 and ferromagnetic spin-1 Bose gases. 
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Introduction. — An understanding of quantum mag- 
netism is important for a large number of phenomena in 
physics. Three well-know examples are high-temperature 
superconductivity, quantum phase transitions and the 
quantum Hall effect. Moreover, it appears that magnetic 
properties will also be very important in another area, 
namely Bose-Einstein condensation in trapped atomic 
gases. The latter has come about because of two in- 
dependent experimental developments. The first devel- 
opment is the realization of an optical trap for atoms, 
whose operation no longer requires the gas to be doubly 
spin-polarized Jl|,||. The second development is the cre- 
ation of a two-component Bose condensate || , which by 
means of rf-fields can be manipulated so as to make the 
two components essentially equivalent [Q] . As a result the 
behavior of both spin-f and spin-1/2 Bose gases can now 
be experimentally explored in detail. 

Theoretically, the ground-state structure of these gases 
has recently been worked out by a number of authors 
and also the first studies of the line-like vortex ex- 
citations have appeared |^,^]. However, an immediate 
question that comes to mind is whether the spin degrees 
of freedom allow for other topological excitations that do 
not have an analogy in the case of a single component 
Bose condensate. It is the main purpose of this Letter to 
show that the answer to this question is in general affir- 
mative. In particular, we show that spinor Bose-Einstein 
condensates have so-called skyrmion excitations, which 
are topological nontrivial point-like spin textures p0[ . 
Having done so, we then turn to the investigation of the 
precise texture and the dynamics of these skyrmions. 

Topological considerations. — To find all possible 
topological excitations of a spinor condensate, we need to 
know the full symmetry of the macroscopic wavefunction 
^(x) = ^/n(x)£(x), where n(x) is the total density of 
the gas, £(x) is a normalized spinor that determines the 
average local spin by means of (S)(x) = C*(x)S ab C&(x), 
and S are the usual spin matrices obeying the commuta- 
tion relations [S a , Sp] = ie a fp SL. Note that here, and in 
the following, summation over repeated indices is always 
implicitly implied. From the work of Ho Q we know 
that in the case of spin-1 bosons we have to consider two 
possibilities, since the effective interaction between two 
spins can be either antiferromagnetic or ferromagnetic. 



In the antiferromagnetic case the ground-state energy is 
minimized for (S)(x) = 0, which implies that the pa- 
rameter space for the spinor £(x) is only [7(1) x S 2 be- 
cause we are free to chooce both its overall phase and 
the spin quantization axis. In the ferromagnetic case the 
energy is minimized for |(S)(x)| = 1 and the parameter 
space corresponds to the full rotation group SO(3). Us- 
ing the same arguments, we find that for spin-1/2 bosons 
the order parameter space of the ground state is always 
equivalent to SU(2) fjj. 

What do these results tell us about the possible topo- 
logical excitations? For line-like defects or vortices, we 
can assume £(x) to be independent of one direction and 
the spinor represents a mapping from a two-dimensional 
plane into the order parameter space. If the vortex is 
singular this will be visible on the boundary of the two- 
dimensional plane and we need to investigate the prop- 
erties of a continuous mapping from a circle S 1 into the 
order parameter space G, i.e., from the first homotopy 
group 7ri(G). Since tt 1 (SU(2)) = wi(SO(3)) = Z 2 and 
7Ti(C/(l) x S 2 ) = Z, we conclude that a spin-1/2 and a 
ferromagnetic spin-1 condensate can have only vortices 
with a winding number equal to 1, whereas an antiferro- 
magnetic spin-1 condensate can have vortices with wind- 
ing numbers that are an arbitrary integer. Physically, 
this means that by traversing the boundary of the plane, 
the spinor can wind around the order parameter at most 
once or an arbitrary number of times, respectively. This 
conclusion is identical to the one obtained prevously by 
Ho [|. 

If the vortex is nonsingular, however, the spinor £(x) 
will be identical everywhere on the boundary of the two- 
dimensional plane and it effectively represents a mapping 
from the surface of a three-dimensional sphere S 2 into 
the order parameter space. We then need to consider the 
second homotopy group 772(G). For this we have that 
ir 2 (SU(2)) = 7r 2 (5*0(3)) = and n 2 (U(l) x S 2 ) = Z. 
Hence nonsingular or coreless vortices are only possible 
for a spin-1 condensate with antiferromagnetic interac- 
tions. It therefore appears that the nonsingular spin tex- 
ture discussed in Ref. ||, is topologically unstable and 
can be continously deformed into the ground state by 
'local surgery' ]l2|j . 

We are now in a position to discuss point-like de- 
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fects. Since the boundary of a three-dimensional gas 
is also the surface of a three-dimensional sphere, sin- 
gular point-like defects are also determined by the sec- 
ond homotopy group 7^(6?) . Such topological excita- 
tions thus only exists in the case of a spin-1 Bose gas 
with antiferromagnetic interactions. We call these ex- 
citations singular skyrmions, although in view of the 
work of 't Hooft and Polyakov it would be justified to 
call them monopoles p3| . For nonsingular point-like 
defects the spinor £(x) will again be identical on the 
boundary of the three-dimensional gas. As a result, 
the configurations space is compactified to the surface 
of a four-dimensional sphere S 3 and we need to deter- 
mine the third homotopy group 773(G). For this we find 
7r 3 (S77(2)) = 7r 3 (50(3)) = K 3 (U(l) x S* 2 ) = Z. Hence 
nonsingular skyrmion excitations exists in all three cases. 

Skyrmion texture. — In view of the fact that magnetic 
gradients are needed to stabilize the spinor condensate in 
the case of a spin-1 Bose gas with antiferromagnetic in- 
teractions J8|, we from now on consider only spin- 1/2 and 
ferromagnetic spin-1 Bose gases, which have only nonsin- 
gular skyrmions as we have seen. Due to the absense of a 
core, the fluctuations in the density Sn(x) = n(x) — (n(x)) 
will be small compared to the average density (^(x)) and 
the energy of the skyrmion can be determined by 

E[Q = ^jdxjdx' 5n(x)x _1 (x,x')<5n(x') 

+ |dx n(x) (|^|VC(x)| 2 - /uB • (S)(x)) . (1) 

Here, x(x, x') is the static density-density response func- 
tion defined by 

£(- V '((^ V ) +167ra ) X(X ' X,) 

= «5(x-x'), (2) 

a is the appropriate scattering length, B is either a ficti- 
cious magnetic field, caused by resonant rf-fields, or a real 
homogeneous magnetic field and /1 is the corresponding 
magnetic moment of the atoms in the trap. Moreover, 
the spinor £(x) can now be conveniently parametrized as 

C(x)-cxp|-^ a (x)5 Q |c Z , (3) 

where S denotes the spin of the atoms and ( z is a con- 
stant spinor that minimizes the Zeeman energy, i.e., it 
obeys £o = S a s if we use the direction of the magnetic 
field as the quantization axis. 

This last equation explicitly shows that the topology 
of the order parameter space is a sphere of radius tt with 
opposite points on the surface identified. If we now as- 
sume that a maximally symmetric shape of the skyrmion 
spin texture is allowed, we can take f2(x) = xu)(x)/x = 



xw(x), with ui(x) a monotonically decreasing function 
obeying w(0) = 2ir and linx c _ s . 0O ui(x) = 0. For this ansatz 
we see that by traversing the whole configuration space, 
we exactly cover the order parameter space once. The 
above spin texture therefore corresponds to a skyrmion 
with a topological winding number 

W =^J dn = ^J d^ k e aPl d^ a d^ 3 d k W (4) 

equal to 1 and located at the center of the trap. 

To obtain the precise spin texture of the skyrmion we 
have to find the function u>(x) that minimizes the en- 
ergy. In general this leads, due to the nonlocality of 
the density-density response function, to a complicated 
nonlinear integro-differential equation, which can only be 
solved numerically. However, we can gain a lot of phys- 
ical insight in the size and stability of the skyrmion by 
the following analysis. Considering first the ideal case of 
zero magnetic field B and solving for the density fluctu- 
ations induced by the spin texture, we have to minimize 
the gradient energy 

E^ d {(} = / dx(n(x))|l|VC(x)| 2 
J 2m 

I dx J dx ' |VC(x)| 2 x(x,x')|VC(x')| 2 • (5) 

For a skyrmion of size A in the center of the trap this 
energy can be estimated to be given by 

•^ = *T(>-(i|Tl) ' ,6 > 

with £ = (167r|a|(n(0))) _1//2 the typical coherence length 
of the condensate. It has a minimum when A/£ is equal to 
((\/33 — 5)/2) 1 / 2 ~ 0.61. Therefore, our maximally sym- 
metric ansatz is indeed justified, even for an anisotropic 
trap, as long as the coherence length is much smaller 
than the size of the condensate, i.e., we are in the so- 
called Thomas-Fermi limit. 

For this result we have assumed the scattering length 
a to be positive. In that case the skyrmion is thus sta- 
bilized by the fact that gradients in the spin texture on 
the one hand lead to an increase in the average kinetic 
energy of the spinor condensate, but one the other hand 
also lead to a reduction of the density, which for repul- 
sive interactions results in a decrease in the average in- 
teraction energy. For a negative scattering length, gra- 
dients in the spin texture lead to an enhancement of the 
density, but this, due to the attractive interactions, also 
reduces the average interaction energy. Therefore, we ex- 
pect skyrmions also to be stabilized in that case. Indeed, 
using the same estimates as before, we find that E siad (X) 
is now minimized for A/£ equal to ((v / 33 + 5)/2) 1 / 2 ~ 2.3. 
However, due to the intrinsic instability of a condensate 
with attractive interactions JL4] , this implies that the size 
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of the skyrmion is essentially equal to the size of the con- 
densate itself and that its spin texture will be strongly 
affected by finite-size effects. In particular, it in general 
becomes anisotropic. 

Note finally, that the above mechanism for the stability 
of the skyrmion excitations is quite different from that 
in the quantum Hall effect near filling fraction v = 1. 
There skyrmions are electrically charged and their size 
is determined by a competition between the Coulomb 
and Zeeman interactions |l5| ]. In the case of a spinor 
Bose-Einstein condensate, the Zeeman interaction plays 
a much less important role and for small magnetic fields 
B <C 87r|a|?i 2 (n(0))/TO^ only directs the spins at large 
distances of the skyrmion. Larger magnetic fields will 
tend to decrease the size of the skyrmion so as to reduce 
the cost in Zeeman energy. 

Skyrmion dynamics. — The most important dynami- 
cal variable of the skyrmion arises from the fact that the 
Euler-Lagrange equations for the skyrmion spin texture 
is invariant under a space independent rotation of the 
average local spin (S)(x) around the magnetic field di- 
rection B [Il6j. Mathematically, this means that if C sk ( x ) 
describes a skyrmion, then exp{— iz9B-S}£ sk (x) describes 
also a skyrmion with the same winding number and en- 
ergy. The dynamics of the variable #(£) associated with 
this symmetry is determined by the full action for the 
spin texture S[(] = J dt(T[(] — E[Q), which contains the 
time-derivative term 



4n(x,f)(*(x,i)ft-((x,t) . 



(7) 



Hence, subsituting ((x,t) = exp{-ii?(i)B • S}C sk (x) 
and making use of the conservation of total particle num- 
ber N to introduce the change of the average local spin 
projection on the magnetic field (A5||)(x) = B-(S)(x)— S 
induced by the skyrmion, we obtain apart from an unim- 
portant boundary term that 



T[Q=n^ J dxn(x,t)(A5||)(x) 



(8) 



Taking now again the density fluctuations into account, 
we finally find that the dynamics of the rotation angle 
is determined by the action 

where (AS 1 ^ *) is the change of the total spin along the 
magnetic field direction and the 'moment of inertia' of 
the skyrmion equals 



I = h 2 J dxj dx' (AS||)(x) X (x,x')<AS||)(x') . (10) 

The importance of this result is twofold. First, from 
the action in Eq. (ft) we see that at the quantum level 



the hamiltonian for the dynamics of the wave function 
t) becomes 



H 



21 lift? 



(11) 



Therefore, the ground state wave function is given by 



*o(0) 



iKt) 



I v 27T, with K an integer that is as close as 



e 

possible to (AS^ *). In this way we thus recover the fact 
that according to quantum mechanics the total number of 
spin-flips associated with the skyrmion texture must be 
an integer. More precisely, we have actually shown that 
the many-body wave function describing the skyrmion 
is an eigenstate of the operator B • S tot with eigenvalue 
NS — K. Note that physically this is equivalent to the 
way in which 'diffusion' of the overall phase of a Bose- 
Einstein condensate leads to the conseration of particle 
number 

Furthermore, the existence of this internal degree of 
freedom becomes especially important when we deal with 
more than one skyrmion in the condensate. In that case 
every skyrmion can have its own orientation and we ex- 
pect the interaction between two skyrmions to have a 
Josephson-like contribution proportional to cos(z9i — #2). 
As a result the phase diagram of a gas of skyrmions can 
become extremely rich and contain both quantum as well 
as classical, i.e., nonzero temperature, phase transitions 
p6fl . In this context, it is interesting to mention two 
important differences with the situation in the quantum 
Hall effect. First, the fact that the spin projection K of 
the skyrmion is an integer shows that these excitations 
have an integer spin and are therefore bosons p9[ . In 
the quantum Hall case the skyrmions are fermions with 
half-integer spin, due to the presence of a topological 
term in the action S[Q for the spin texture Second, 
in a spinor Bose-Einstein condensate the skyrmions are 
not pinned by disorder and are in principle free to move. 
Both differences will clearly have important consequences 
for the many-body physics of a skyrmion gas. 

Focussing again on the single skyrmion dynamics, we 
can make the last remark more quantitative by using the 
ansatz <^(x, t) = £ sk (x — u(t)) for the texture of a mov- 
ing skyrmion, which is expected to be accurate for small 
velocities du(t)/dt. Considering for illustrative purposes 
also only the case of a maximally symmetric skyrmion 
near the center of the trap with u/£ <C 1, we find in a 
similar way as before that the action for the center-of- 
mass motion of the skyrmion becomes 



S[u] 



dt -M 
2 



du(t) 

dt 



(12) 



where the mass, which in the more general anisotropic 
case is of course a tensor, is now simply given by 



W = ^- / dx / dx' X (x,x')(v s )(x) ■ (v s )(x') (13) 
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in terms of the superfluid velocity of the spinor conden- 
sate (v s )(x) = —ih( sk (x)V( sk (x)/m. The skyrmions 
thus indeed behave in this respect as ordinary particles. 
It should be noted that in contrast to Eq. (^) there ap- 
pears no term linear in du(t)/dt in the action S[u]. This 
is a result of the fact that we have performed all our 
calculations at zero temperature. In the presence of a 
normal component, we anticipate the appearance of such 
a linear term with an imaginary coefficient. This will lead 
to damping of the center-of-mass motion of the skyrmion. 
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